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Abstract. We rigorously derive multi-pulse interaction laws for the semi-strong interactions in a 
family of singularly-perturbed and weakly-damped reaction-diffusion systems in one space dimen- 
sion. Most significantly, we show the existence of a manifold of quasi-steady W^-pulse solutions 
and identify a "normal-hyperbolicity" condition which balances the asymptotic weakness of the 
linear damping against the algebraic evolution rate of the multi-pulses. Our main result is the 
adiabatic stability of the manifolds subject to this normal hyperbolicity condition. More specif- 
ically, the spectrum of the linearization about a fixed A^-pulse configuration contains essential 
spectrum that is asymptotically close to the origin as well as semi-strong eigenvalues which move 
at leading order as the pulse positions evolve. We characterize the semi-strong eigenvalues in 
terms of the spectrum of an explicit N X N matrix, and rigorously bound the error between the 
Af-pulse manifold and the evolution of the full system, in a polynomially weighted space, so long 
as the semi-strong spectrum remains strictly in the left-half complex plane, and the essential 
spectrum is not too close to the origin. 



1. Introduction 

It is not uncommon that evolutionary partial differential equations have finite dimensional sub- 
manifolds which are approximately invariant, and robustly stable, in the sense that initial data 
which starts close to the sub-manifold remains close. For these proximal orbits, the full dynamics 
typically can be reduced to a finite dimensional system, often with nontrivial dynamics of its own. 
Asymptotic analysis is adept at constructing such manifolds; determining the stability of the sub- 
manifold is more problematic, as the natural object which arises in the linear stability theory 
may have non-trivial time-dependence. The road from a time-dependent linear operator to the 
properties of its semi-group is long and hard. A tractable sub-case arises when the flow on the sub- 
manifold is slow in comparison to the exponential rates which characterize the decay of proximal 
orbits towards the sub-manifold; in this case we say the sub-manifold is normally hyperbolic. It 
is natural to investigate the case in which the normal hyperbolicity is lost due to a decrease in 
the exponential decay rates of the sub-manifold. We address this question within the context of 
semi-strong multi-pulse interactions in a weakly-damped reaction diffusion system. 

The study of pulse interaction in diffusive systems has a long history, in particular the activator- 



inhibitor systems modeled by the Gierer-Meinhardt 27 and the Gray-Scott 31 equations have 
spawned a substantial literature. These reaction-diffusion systems are comprised of two chemicals 
which feed an autocatalytic reaction that drives pattern formation. Our analysis is particularly 
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motivated by the Gray-Scott Model 



(1.1) 



Ut + A{1 -U)- UV^ 

Vt^DV,.,-BV + UV^, 



where U and V denote the concentrations of the chemical species. The semi-strong regime of the 
Gray-Scott equations, presented in 12 13 15], occurs when the U component experiences strong 



diffusion, but is weakly damped and strongly forced. In particular, there is a (nontrivial) balance 
which maintains the two components at an 0(1) level within the spatialdomain where the slowly 
diffusing component is present. This leads to a rich stability structure in which the localized 
species manifests a long-range interaction through the delocalized species. This simple system 
affords an ideal arena for the study of reductions of infinite dimensional dynamical systems to 
finite dimensional sub-systems. This paper studies a generalization of the rescaled system. This 
generalization maintains the aforementioned balance while enjoying a variable rate of linear damping 
and a homogeneous polynomial nonlinearity in U = {U,V)'^ , 



(1.2) 



Ut = 
Vt = 



K.X - V 



=:T{U), 



where an, a2i > and /i > 0, a > 0, p > 0, ai2 > 1 and Q!22 > 1- For < e ^ 1 the system affords 
a natural competition between long-range, U, and short-range, V, interactions. Such competitions 
are wide-spread in physical settings, arising for example from the balance between entropic and 
electrostatic interactions in ionic solutions which drive morphology generation in solvated charged 
polymers, 
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We address the existence, dynamics, and particularly the adiabatic stability of A^-pulses in the 
semi-strong interaction regime. This regime features pulse-like structures in the localized species 
whose positions, amplitudes, and stability all evolve at leading order subject to an effective mean- 
field generated by the pulses coupling to the delocalized, rapidly-diffusing species. The novelty of 
our study lies in the asymptotic approach of the essential spectrum to the origin, characterized by 
the weak linear-damping rate, e"^, which may lead to a loss of normal hyperbolicity (adiabaticity) 
as the slow decay associated to the essential spectrum competes with the perturbations generated 
by the pulse evolution. 

For fixed TV G N+, we rigorously derive the existence and adiabatic stability of semi-strong 
A^-pulse configurations. This study requires a careful analysis of the linearization about the A^- 
pulse configurations, including resolvent and semi-group estimates as well as a characterization of 
the point spectrum. Generalizing prior results, 12 14 36 , we show that there is a set of point 



spectra, the semi-strong spectrum, which evolve at leading order in conjunction with the localized 
pulse configuration, and can be characterized in terms of the eigenvalues of an explicit N x N 
matrix. Extending these results, we rigorously show that so long as the finite-rank spectrum remains 
uniformly within the left-half plane and the linear damping is stronger than a critical value, then 
the iV-pulses generate an adiabatically stable manifold which affords a leading order description of 
the pulse interactions. A key element of the analysis is the development of semi-group estimates, via 
a renormalization group (RG) approach, corresponding to the evolving pulse configurations. This 
result requires a form of normal hyperbolicity for the manifold of semi-strong iV-pulses, balancing 
the flow on the manifold against the weak linear decay, see Theorem |l.l| and following discussion. 

There is a developed literature on the stability of viscous shocks and traveling waves, for which 
the essential spectrum touches the origin, see 24 35 43 44 53 54] and references therein. We 
emphasize two important distinctions between our results and these two bodies of work. For the 
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nonlinear conservation laws, |35[|53[|54| , the stability estimates do not "close" in the following sense: 
initial perturbations are considered in a space, L^(M) for example, which is not controlled by the 
estimates at later times - the decay is in a norm which does not control the initial perturbation. 
As a consequence the process cannot be iterated; one cannot restart the perturbation analysis at a 
later time in the flow. For the traveling waves under essential bifurcation, 24 43 44 , the solution 



oscillates temporally in a neighborhood of a fixed structure: the system is controlled by a single, 
temporally fixed linearized operator, thereby avoiding the issue of competition between decay rates 
and secular forcing arising from a time-dependent linearization. Neither family of results extends 
trivially to encompass the semi-strong interactions, in which the underlying structure, the TV-pulse 
positions and amplitudes, evolve at leading order generating concomitant changes in the associated 
linearized operators. In the setting of ( 1.2 1, one is faced either with developing semi-group estimates 
on a time-dependent family of operators, as in 45 , or in the approach we follow here: taking the 



linearized operators to be piece-wise constant in time with a renormalization of the flow at each 
jump in the linearization. This later, iterative approach requires estimates on the decay in the same 
norm used to control the initial perturbations. 



1.1. Prior results on semi-strong interactions. The semi-strong interaction regime is interme- 
diate between the weak interaction regime and the strong interaction regime. In the weak regime, 
the pulses are localized in each component, which returns to its equilibrium value between adjacent 
pulses. The mutual interaction of localized structures is exponentially weak in the pulse separation 
distance, and consequently there is no leading order influence of pulse location on the shape or 
the stability of the pulses. The weak interaction regime has been well-studied in reaction-diffusion 
systems, see 19 20 41 42 . The strong interaction regime arises when the pulses are sufficiently 



proximal that the values of their localized components do not return to equilibrium values between 
the pulses. This leads to self-replication, collision, annihilation, and other strongly non-adiabatic 
behaviors. There has been little theoretical investigation of the strong interaction regime, which 
is typically investigated using numerical techniques. In the semi-strong regime the pulses have 
both localized and delocalized components, with the delocalized components varying slowly over 
the support of the localized components. This regime has been studied both formally 12 16 36 



and rigorously 17 47 



In 36 



These previous works have largely focused on the Gierer-Meinhardt and the Gray-Scott models, 
a formal study of the TV-pulse semi-strong interaction regime for the generalized Gierer- 
Meinhardt model was presented. In particular, expressions for the semi-strong spectrum and the 
ordinary differential equations for the dynamics of A^-pulse configurations were derived. In 16 , 
a general system that includes both the Gierer-Meinhardt and the Gray-Scott model was studied. 
The semi-strong two-pulse interaction was investigated, where formal results for the asymptotic 
stability were determined in particular regimes, along with ordinary differential equations governing 
the dynamics of pulse positions. The extension of these results to the A''-pulse was also discussed. 
However, the derivation of conditions under which the iV-pulse manifold is adiabatically attractive 
under the full fiow of the PDE is outside the scope of the prior work. 

In 17|, the 2-pulse semi-strong interaction regime was rigorously studied for the regularized 
Gierer-Meinhardt model, which has a strong linear damping rate, extending the renormalization 
group approach developed in 41 to obtain appropriate semi-group estimates on families of weakly 
time-dependent linear operators. In 47 , the renormalization group approach was used to establish 
the adiabatic stability of 2, 3, and 4-pulse configurations within the semi-strong interaction regime 
for a three-component system with two inhibitor components and one strongly-linearly-damped 
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activator component. The renormalization group techniques have also been used to study quasi- 
steady manifolds in noisy systems, [32|, and coupled dispersive-diffusive systems, 40 . 



The homogeneous nonlinearity considered in (1.2) is characterized by the quantity 
(1.3) 



:= ail 



ai2a2i 



a22 - 1 

In particular, the admissibility conditions typically require 9 < 0, see Remark |1.1| The scaled 
Gray-Scott system, in 



12 



13 



exceeds the critical value, a 



, corresponds to 9 = —1, however its linear damping rate, a = |, 
permitted by Theorem ] l.l[ This suggests that the evolution of its 
iV-pulse solutions may not be describable by an A'^-dimensional system parameterized solely by the 
pulse locations. It may be necessary to include the impact of the essential spectrum - which may 
manifest itself as a long, low shelf (asymptotically small in L°° but large in L^) which surrounds 
the pulse region. A detailed description of this supercritical regime is an intriguing open problem 
that is outside the scope of this work. 




JX- 

PJ 



Figure 1. A typical iV-pulse configuration for the coupled system (1.2) 
V component is localized at the pulse positions pj 



The 



The U component has an 
approximately constant value, g^, on the narrow pulse intervals, and is slowly 
varying in between the pulses, reaching its equilibrium value of 0(e^"/^) as a; — >■ 
±oo. 



1.2. Presentation of the main results. Our first result is the existence of the semi-strong A^- 
pulse configurations parameterized by the localized pulse positions, p = {pi, ■ ■ ■ ,pn)^ G M^. In 
general, a manifold formed by the graph of a function $ = $(p ) is invariant under a flow Ut = J'iU) 
on a Banach space X if and only if 

(1.4) (/-^p~)^(<f(p))-0, 

where tt^ is the projection onto the tangent plane of the manifold at p. Indeed, the classic manifold 
formed by the translates of a traveling wave solution is invariant precisely because its wave-shape 
satisfies the so-called traveling wave ODE, 

= J-($), 

for wave speed c. The manifold is parameterized by the translates A4 :— {$(• ~ p)\p G K} of $, 
and the projection off of the tangent plane of A4 has in its kernel. In cases not arising from a 
natural symmetry of the flow, such as translation or rotation, construction of an invariant manifold 
is a non-trivial endeavor, see (5^ for example. 
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We do not perform this calculation, rather we construct a manifold with boundary, which is 
approximately invariant. More specifically, we construct a manifold for which the left-hand side of 



(1.4) is sufficiently small in an appropriate norm, and which has a thin, forward invariant neigh- 



borhood which attracts the flow from a thicker neighborhood - at least up to times that the flow 
hits the boundary of the manifold. In Theorem |2.1[ we reduce this construction to the solution of 
an N X N system of nonlinear equations which connects the positions of the localized pulses, p, to 



the amplitudes of the delocalized component at the pulse position qi = U{pi) for i 
family of solutions, q = q{p) gives rise to the semi-strong iV-pulse configuration 

^iix;p] 
<S>2{x;p] 



1, 



,N. The 



(1.5) 



defined in (2.9|. In Proposition 4.2 we characterize the spectrum of the linearization Lg of (1.2) 



about $, showing that the point spectrum consists of N eigenvalues localized near the origin 
together with the semi-strong spectrum, denoted crss(p), which can move at leading order as the 
pulse positions p evolve. Moreover, the locations of the semi-strong spectrum are determined by 



the eigenvalues of an x iV matrix Afx defined in (4.20). This motivates the following definition 



c 
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Figure 2. The contour C introduced in (1.6). The essential spectrum of the 
linearization Lp about an iV-pulse solution, see (4.1), extends from (— c», — ^e"] 
(orange line) while the point spectrum consists of N broken translational eigen- 
values near the origin (orange X), with associated spectral projection n^, together 
with a collection of semi-strong spectrum (orange o) which move at leading or- 
der as the pulse configuration evolves. For an admissible pulse configuration, the 
semi-strong spectrum lies to the left of the contour C. 



Definition 1.1. Fix e > 0. An open, connected set K. C ffi is an admissible family of N -pulses if 
there exist b,C,£ > and v € (0, [i) such that the following hold. 

(1) For all p ^ JC the semi-strong spectrum crss(p) lies to the left of the contour C C C, of the 
form 

(1.6) C=CUC_UC+, 

where Cy = |— e"i/ -I- is|s e [—6, 6]} and C± = |— e"z/ ±ib -\- se^^ \s £ [— oo, 0]|, see Fig- 
ur^^ 
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(2) For all \ eC and p E K, 



(1.7) 



< C 1 



e 



where k\ e\J A + e"/i is a scaled distance of X to the branchpoint of the essential spectrum 
of Cp. 
(3) K- C K-i where 



Pi+i 



p, >£|lne|} 



(1.8) /C^:={peM^ 
fori = l,--- ,N-1. 

To each admissible family of pulse configurations JC C M^, we associate the A^-pulse manifold 

(1.9) M := {$(-;p) pe /c} . 

Our first result is that a non-trivial portion of the semi-strong A^-pulse interaction regime is admis- 
sible if the corresponding single-pulse is linearly stable. 



Proposition 1.1. Admissibility. Consider the system (1.2). Let the associated single-pulse 



solution defined by (2.9) be linearly stable, that is, if except for a simple translational eigenvalue at 



the origin, the point spectrum of the linearization about a single-pulse lies to the left of a contour 
Cu.b for some i>,b > 0. Then for every N e N+ there exists Eq > and such that for all S > 
sufficiently large and for all < e < the portion of the semi-strong domain satisfying 



(1.10) 

is admissible. 



IC 



\pe 



Api > Je-(i+"/2) for i = 1, 



,iV- 1 



Remark 1.1. The single-pulse can be linearly stable only if 9 < 0. For 9 > the linearization about 
the single pulse has a real eigenvalue which lies between the origin and the ground-state eigenvalue 
Xq > of the operator Lq defined in (4.23). 
which are admissible, even when the underlying one-pulse 



However there are semi-strong two-pulse regimes 
not linearly stable, see 12 . The 



IS 



restriction on S in (1.10 1 specifically precludes the oscillatory instabilities which can arise through 
pulse-pulse interactions, see 13 and Theorem 1 and Figure 3.1 of 17 for examples of semi-strong 
spectra crossing into the right-half complex plane as the pulses become too close. Pulses of this type 
typically have complex spectra; an examination of real eigenvalues is not sufficient for stability. 

The pulse dynamics are driven by the projection of the residual -^('I') onto the iV-dimensional 
space associated to the broken translational eigenvalues clustered within 0(e'') of the origin, where 

measures the strength of the tail-tail interaction of the localized pulse component, and the value of 
r can be taken arbitrarily large by increasing £ in ICi. The spectral projection onto the translational 



eigenvalues, irp (see (4.37)), induces the complementary spectral space 



(1.11) 



Xp - 



{u\\\u\\ 



X <00, TTpU 







where the X-norm, defined in (1.25), is locally polynomially weighted about the position of each 
pulse. Our main result is the adiabatic stability of the admissible semi-strong iV-pulse manifolds 
in this norm, so long as the asymptotic damping is not too weak. 
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Theorem 1.1. Adiabatic stability. Consider the system with 6^0 satisfying the normal 

hyperbolicity condition a £ [0, ^). Let JC be an admissible domain of pulse configurations and fix 
v S [0,/i). Then the associated N-pulse manifold Ai is adiabatically stable in the norm \\ ■ \\x, 
up to 0(e^~"). That is there exists Mo,cJo,eo > such that for all e < Cq and any initial data 
Uo ■■= iUo,VoV of the form 

(1.12) Uo^<i>ix;p) + Woix), 



with dist(p,9/C) > do o.iT'd \\Wo\\ < Moe"|lne| ^, the solution of (1.2) can be uniquely decomposed 
for t G [0, Tf,] as 

(1.13) U{x,t)^<S>{x;p{t)) + W{x,t), 

where p £ K, is a smooth function of t, the remainder W{-,t) e Xj;(^t) satisfies 

(1.14) \\W{t)\\x < M (e-'"-'\\Wo\\x + e^-") , 

and the time Tf, for the pulse configuration p — pit) to hit dIC satisfies Tf, > ModQC^^. Moreover, 
during this time interval the pulse configuration evolves at leading order according to 

(1.15) ^ eQ-U{p)q' + 0{e'-'^,e\\W{t)\\x,\\W{t)\\j,) , 

where Q{p) = diag((f(p)) is the diagonal matrix of amplitudes, and the anti- symmetric matrix A 
is given by 

I + q/2 



"21 00 



(1.16) [A],,= ""^ { k = j, 



in terms of p and the pulse profile (f)o of (2.4) through the total mass of its powers, see (1.291 



Remark 1.2. Since the matrix A has negative entries above the diagonal, and positive below, 
the pulses generically repel each other, with the right- and left-most pulses move right and left 
respectively. Moreover, if the pulse configuration enters the weak regime, then then the inter-pulse 
separations are monotonically increasing and the time Tb is infinite. 

Remark 1.3. The error from the \\W\\x term may initially dominate the pulse motion if \\Wa\\x 
is initially 0(e"). However, after at^i 0{e~°') transient, during which the pulses move a negligible 
0{e°') distance, this term becomes higher order. After this transient the error bound on the pulse 



evolution ( 1.15[) is sharp when the pulses are uniformly in the semi-strong regime, that is each pulse 



separation satisfies age^"''"/^ < Api < oie^"*""/^. However, if one or more pulse separations become 
large, and some pulses enter into a weak regime, then the leading order pulse evolution is still given 



by the first term on the right-hand side of (1.15). A proof of this requires the decomposition (6.1) 



to be built upon a linearization about the more accurate pulse ansatz, $ + <i>c . Thus we must analyze 



the linearization about $ to construct $c; see (6.5), and then linearize about $ + <i>c to perform the 



RG iteration. For brevity of presentation, we have forgone these technicalities. 

The paper is organized as follows. In section 2 we present the existence of the semi-strong A'^- 
pulses, whose graphs form the semi-strong A'^-pulse manifolds. In particular, we uniformly bound 
the components of q both from above, and away from zero from below. In section 3 we develop 
estimates on the resolvent of a key linear operator and on the residual, obtained by evaluating the 



right-hand side of (1.2) at a semi-strong A^-pulse $. In section 4 we characterize the spectrum 
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of the linearization L^'' about an admissible A^-pulse, and in section 5 we develop resolvent and 
semi-group estimates on the full operator. In section 6 we use the renormalization group approach 
to develop nonlinear estimates on the full system and obtain the adiabatic stability results. Section 
7 presents technical estimates used in section 6. We conclude with a discussion which motivates a 
more general relation between the normal hyperbolicity and adiabatic stability. 

1.3. Notation. We fix the number N e N+ of localized pulses and the minimal pulse separation 



parameter, £ > 0, defined in ( 1.8 1. For each fixed pulse position p G /C we define the following norm, 
(1-17) WfWwl^^Ufh^ + Wd.fh^. 

where ^ is a smooth, positive, mass one function with support within (— ^/4,^/4). We also define 
its p translates 

(1.18) ^y.^^{x-p,). 
The W^'' norm controls L°° , since for any x,y G M, 

(1.19) \uix)\<\uiy)\+ Jju'{z)\dz, 

and multiplying this inequality by the mass-one function ^(y), and integrating over y yields 

(1.20) \u{x)\ < Wu'Wli + £ \ayHy)\dy = hllv^^i.i. 

For each /3 > and pulse configuration p e /C wc introduce the locally weighted space defined 
through the partition of unity 

(1-21) {xAp)}^=i, 

which is subordinate to the cover {{sj^i + l,Sj — l)|j = 1,...A^} with sj = ^h±^zPz^ for j = 
{l,...,iV — 1} while So — — oo and sjv = oo, see FigurejSj The space is defined through the 
corresponding norm 



N 

(1-22) \\fh.^^^ = J2\\{l + \x-p,f)x,f\\,,, 

which controls long-wavelength terms, uniformly in each Xj window. In particular, we have the 
estimate 



(1-23) Xjf < C\\f\ 



where the hat denotes the Fourier transform and W^'°° is the classical Sobolev norm. We define 
the windowing {fj ~ Xjf}jLi of ^ function / with respect to the pulse configuration p. This affords 
the decomposition 

N N 

(1-24) / = E/^ = E/^^- 

The main results are stated over the Banach space X defined by the norm 

(1-25) II^^IU-||/i||wi,i + ||/2||hi., 



Adiabatic Stability under Semi-Strong Interactions 



9 







Xj-1 








U 


I 






\ \ 




-1 






j 




^1 



Figure 3. The localized component of the pulse solution associated partition of 
unity {xj}- 



for F = (/i, /2)^ where 7 = maxlan, 021}, and 

(1-26) \\f\\i^^--=\\r\\h + \\f\ 



2 



The H}^ I norm controls the usual norm and affords the nonlinear estimate (3.21) required to 



control the nonlinearity X, as in (7.8), for the case 0:12 = 2. 

For G G L2(i^)fcx! ^-^j^ components [G],j = gij G (K) , we denote the tensor operator (g) G : 
L^(R) I— >• M'^^', which acts on h & L^(M) by component- wise inner product, 
(1.27) [^G-hl^^^{g,,,h)^,. 

If F e L'^{M.y^'^ then the tensor product F ig) G is a finite rank map that takes h e L^(M) to 
F ^ G ■ h G L'^iM.y^™ through the usual matrix multiplication of F with ®G-h. In particular, for 
each p € IC we define the associated windowing tensor 

(1-28) ® = ® ixi,--- ,Xn), 

where {xj}jLi is the partition of unity associated to p and the superscript T denotes transposition. 
The windowing tensor reproduces the masses of the windowings of /, in vector form, 

where we denote the mass of ft, € L^(IR) by 

(1.29) h:= h{x)dx. 

Jr 

Combining the windowing tensor with the vector function 

(1-30) r = i^i.--- ,^n), 

where the {^j jjLi are defined in ( 1.18 ), yields the windowing tensor product £, (E) x which is a rank 
N map from L'^{M.) to L'^{R) defined by 

N N 

(1-31) ® X • / = J2^XjJ)2U^) = ^7, ^x). 

The windowing tensor product replaces / with a sum of N compactly supported functions, centered 
at the pulse positions p with the same mass in each window as /. Indeed, we have the equality 



ix^ ■ (/ - e ® • /) = 0. 
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For vectors q £ M and /3 £ M wc introduce the component-wise exponential, 
(1.32) 



2. Construction of the Semi-strong TV-Pulse Manifolds 



The localized component of the pulse, 0Oi see ( |2.4[ ), decays exponentially at an 0(1) rate. Sub- 
sequently we take the pulse-pulse separation distance £ > in /C^, recall ( |1.8[ ), so large that the 
locahzed tail-tail interaction (j)^{{pj — Pj+i)/2) ~ 0^(^/2) ~ 0{e^) for some r — r{i) ^ \n£ > 2. 

To each vector of amplitudes q G 'R^ we will associate a semi-strong iV-pulse configuration 

<f> — {^i{x;p,q),^2{x;p,q))^, see Figure[l] Moreover we slave the vector q to the positions p of 
the localized pulses via the mean-Geld equations 

(2.1) q, = $i(p,), 

for j = 1, . . . , A'^. The construction begins with the localized fast-pulse the unique solution of 
the equilibrium equation 

(2.2) =0, 

which is homoclinic to the origin and symmetric about x = Pj- The existence of for a22 > 1 is 
immediate as the equation has a first integral. A simple re-scaling shows that (jjj can be written in 
the from 

(2.3) 

where 0o is the unique homoclinic solution to 
(2.4) - 0^' - ^0 

which is even about a; = 0. The second component of $ is the sum of the localized pulses, 



If \ — ^21 /{a22 — 1) / / \ 

<Pj[x;Pj,qj) = qj '(po(x-pj), 



ia.li 



(2.5) 



N 

$2 = ^(/)j(x;p,g). 

The first component of $ is the mean field generated by the localized pulses, obtained by solving 
the first equation of (1.2) at equilibrium with U replaced in the nonlinearity by the local constant 
value qj when acting on (^j. The result is a linear equation for $i, 

N 

(2.6) - e"A**i - Ij'^T" + ^"^'z' = 0' 

Introducing the operator associated to the essential spectrum, 

(2.7) LI, = -e-^dl + eV, 



a/2P 



with inverse denoted L^^^, we rewrite (2.61 as 

(2.8) L^i (^$1 - e 
and define the A'^-pulse configurations 

(2.9) Hx,p,q) 



N 



'^i{x,P,q) 






_ ^2{x,p,q) 







E,=i<?^j(a;) 
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For A e C we introduce the Green's function G\{x) associated to L^i + A which enjoys the 
property 

(2.10) m, + X)-'f={Gx*f){x). 
From the Fourier transformation, we determine its exphcit formula 

(2.11) Gxix) ^. 



'kx\x\ 



where kx was introduced in Definition |l.l[ A central role is played by the scaled two-point correlation 
matrix Gn{^) of the Green's function ( 2.11[ ) 

(2.12) [Qn{p,X)],, 

where p G JC for some minimal pulse separation i > 0. The unsealed version of the two-point 
correlation matrix is denoted 

J2 



(2.13) 

The following theorem shows the existence of solutions to the mean-field equation. 



Gn{j>,>^) = ^^Ar(p, A). 



Theorem 2.1. (Existence of N -pulse configurations) Fix the pulse separation £ > and 
N e N+. There exists a unique function q = q{p) for which ^ (x,p,q{p)) solves the mean-field 
equation (2.1). Moreover, the solution takes the form if — (fo(p) + ^'^^"^^Qiip) , o,nd, is uniformly 
bounded, component-wise from above and away from zero. In particular, qo solves 

(2.14) 

and q admits the expansion, 



^-^gN{p,Q)%' = -i-e"^%. 



1 



(2.15) 



9(p) 



P 



r^jv'(P,0)l 



O 



{'"") ■ 



where the exponential is component-wise with 9 defined in (1.3). 



Proof. An application of (2.101 with A = to the first component of (2.9) evaluated si x = pk 
yields the equation 



(2.16) 



N 



where, for A = 0, the Green's function takes the form 



(2.17) 



Go(x) 



The system of equations in (2.16) may be written in the vector form 
(2.18) 

where R € M.^^^ has entries 

(2.19) R,fe=e"/2-i(Go*0^"-)(pfc). 
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Substituting for Go from (2.17), for from (2.3), and recalling the definition (1.3) of 9, we obtain 



(2.20) 



R 



from which we see that R has an 0(1) limit as e tends to zero. For j ^ k the exponential decay of 
4>Q permits us to Taylor expand the exponential about y = Pj, 

(2.21) e-^'^^v^lP^-^l ^ e~^'^*v^lP^--P^I-(7jfee^+^VAte-''''*^|P'-"P^I(y-Pj) + 0(e2+"|z/-Pjf ), 



where cjjk — sign(pfe ~Pj)- Substituting this into (2.20) and evaluating the leading-order integral, 
we have the expansion 



(2.22) 



R 



'3^ I — 7 



where the ©(e^"*""/^) terms are zero due to parity, and we are applying the total mass notation 
/ := Jjg f{s)ds. For the case j = k we consider the Taylor expansion of the exponential for y < pj 
and for y > pj. The 0{y — pj) terms do not integrate to zero and we obtain 



(2.23) 



R 



Therefore, using both (|2.22|) and (|2.23|), we see that the equation (|2.18|) can be written as 
(2.24) 



ew(p,0)+ei+"/2Ri(p,e) 



where Ri e 



is uniformly bounded for |e| < eo and smooth in p. Rescaling 



removes the e dependence of the exponentials. The existence of the solution to (2.14) of the form 
(2.15) follows from the implicit function theorem and the invertibility of Qn{p, A), indeed 

N-l 



(2.25) 



detg^(p,A) = n (l 



-2fcA(Pi + l-Pi) 



which is non-zero. Moreover, the inverse oi Qn is tri-diagonal and can be constructed explicitly, 
see 36 for details. The positivity and uniform bounds on q for p in the semi-strong regime follow. 
The existence of the solution q to (2.24) then follows from a perturbation off of ifo. □ 



3. Bounds and Residual Estimates on Semi-Strong TV-pulses 

Fix a set K. of semi-strong N pulse solutions as constructed in Theorem |2.1[ with minimal pulse 
separation distance £ > 0. For all p € /C, in subsection 3.1 we establish estimates on (Lfj^ + A)"""^/ in 
various norms, in subsection 3.2 we obtain bounds on the norms of the semi-strong A^-pulses $(p), 
in subsection 3.3 we establish a result allowing for subsequent reduction of a finite rank operator, 
and in subsection 3.4 we establish estimates on the residual F{^{p)). 
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3.1. Linear estimates. We recall Lf^ and kx introduced in (2.7) and Dcfinition |l.l| respectively. 
For each p G K,, we have the partition of unity {xj}jLi given in (1.21 ) and the weighted-windowed 
norm - defined in ( 1.22 ). 

Lemma 3.1. There exists C > such that for any f e Li(M) or f € W^'{R) and for any X G C 
\ (— cx), — e"/i), the following estimates hold: 

WfWwy] 



(3.1) 



\\{Ll,+X)-'f\\^y < 



mill { ll/IUi, 



\kx 



(3.2) ||(L?i-f A)-Vl|L^+|5i(fcA)|||(i?i + A)-VlUi < ^^ll/IU^' 

(3.3) |ia,((L?i + A)-\/)|U~ + |5R(fcA)|||9.((i?i + A)-i/)|Ur < Ce^WfU^, 

where 3? denotes the real part. Moreover, for all p E IC, / G L\ p>(M), and A G C \ (— oo, — e"/i) 
have the small-mass estimates 

(3.4) \\{Ll,+\)-'f\\^.,. < C^{\®x-f\ + \kx\\\f\\LQ, 



(3.5) 



ll(in + A)-Vll 



1 



I®x-/I + II/IIli 



The small-mass estimates are useful in section 6 when we examine the difference of two linear 
operators, where their difference is large, but the difference has small mass in each window. 

Proof. We introduce g{x) := {L\i + A)^^/ = {G\ * f){x) where the Green's function Gx is given in 
(2.11 ). From the identity g' — G'^ * f and convolution estimates 34 , we have the bounds 

^2 



(3.6) 
(3.7) 



IWh^ <\\G'xh4f\\L^ <C^^\\f\\L^. 

Ugh^ < Uh^iGx * /)(x)|U^ < ||G,|U^.|l/|Ur < c£j\\f\\L^, 



for ^ defined in (1.181. These two estimates establish the bound in (3.1). For the W^' bound 
we first observe that 

,2 



(3.8) 



WWl^ = \\{Gx* f'){x)h^ <\\Gxh4f\\L^ <C 



\kxmkx 



-11/11 



w 



Combining this with the estimate 



(3.9) 



I1^5I1li < llCllLi|l(GA*/)(a;)|k^ < IIGaII 



< c 



\kxmkx) 



11/11 



yields (|3Tj). 

To establish (3.2) we observe that 

(3.10) 



llfflU- =|1(Ga*/)(x)|U^ < ||Ga|1l^|1/|1li <G^||/|Ur, 



Ml^ =|l(GA*/)(x)|Ur < |1Ga|1li|1/I1li <G- 



(3-11) nv". ■ ..v-.n.- _ n^.n.-n.n^- - - l^^jSjj^^^ 

The bounds (3.3) follow from similar estimates applied to g' = G'^ * /. 



■11/11 
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For the small-mass estimates, we window / through its partition of unity, as in (1.24 1, so that 
(3.12) g, = Gx * fj 



satisfies g = gj. From the definition (1.22) of the windowed norm we see ^ = X]j II/jIIl} ■ 

We decompose each into a smooth, localized term and a massless part, 

(3.13) /, = /,e,+y;, 

for Uj G W^''^{M.) and defined in (1.18). Clearly for any /, ||/||li < II/IIlJ .■ Next, we examine 

(3.14) = / [dAx-p,)]\y,\dx< I \{x-pj)y]\dx = C f \{x - p,) {f, - f,^,) \dx, 

Jk Jr Jr 

(3.15) <c(^mLl^ + \\f\\L^ Jj{x~-p,)^,\dx^ <C||/,L; . 

We decompose gj = gj^i + gj,o where gj^i — fjG\ * S^j and gj^ ~ G\ * y'^ = G\ * yj. Estimating gj^i 
using (2.10) and ( |3.1[ ), we have 

(3.16) 



The function G'^ has a jump at a; = 0, so that = [G"] + €^60, where [G'l] is the point-wise second 
derivative of Gx and Sq is the delta function at x = 0. This yields the estimate, 

HJl^ < (II [G'i] h^+e') l|y,IU. < ^jII^^'II/.IIl;,. 



(3.17) 



Using (3.3 1, we find 

(3.18) U9,..o\\l^ < lieilLi||G^*z/,|U^. <G||G^|Ui||y,|Ui <Ge2||/,||^;^.. 

Summing over j, we have (3.4). The inequality (3.5) follows using (3.3 1 and (3.2) respectively. □ 

3.2. Bounds on ^(p ). The following lemma establishes bounds on the A^-pulse solutions over each 
admissible set. 



Lemma 3.2. Let K. denote a family of N -pulse configurations, as constructed in Theorem 2.1 
Then for all /32 > there exists a constant G > such that VpG /C and k = 1, ■ ■ ■ 



(3.19) 



(3.20) ||9.$i||l- + ||5p,$i||L, 



l*i| 



-a/2 





$1- 








+ 



-a/2P 



<Ce, 



where x = x(':P) ih,e partition of unity subordinate to p defined in (1.211. Moreover, for all 
/3 > and all u e - there exists G > such that 



(3.21) 



,<Ce\\v\\^. ^. 



Proof. To establish the bounds on the first two terms in ( 3.19| ), we apply (3.2), with A = 0, to (2.8) 
and recall that the q are uniformly 0(1). For the final term of (3.19), we take dp^ of (2.8) to obtain, 



(3.22) 



dqi 
dpk ' 
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Inverting L^i and applying (3.3) with A = to the first term on the right-hand side, and (3.2) 
together with the bound on dp^qj to the second term establishes the result. The first estimate in 
(3.20) follows from applying (3.3) to the derivative of (2.8) at A = 0. The smooth dependence of q 
on p and the change of variables p = pe^+"/^ imply that | = 0(e^+"/^), from which the second 
and third estimates of (3.20) follow. The final estimate of (3.20) follows from the 0(e) L°°-bound 
on dx^i, the fact that ^<I>^^ — x ■ '7*'^^^ (Pj) — 0, and the exponential decay of $2 away from x 
The estimate (3.21) follows from the bound 

Xj 



Pj- 



1 



Po\ 



< Ce. 



□ 



3.3. Finite-rank reduction lemma. The following Lemma is a key tool in identifying the finite- 
rank reduction of the singularly scaled operator Lf^ + A, by describing its structure. 

Lemma 3.3. For each pulse separation parameter £ > there exists C > such that for a^/ A G C 
\ (—00, —e^fi) and p € )Ci, the following holds 

(3.23) {iLl, + \)-'f,g)^,-i®x-ffGM{p,X)^X-9 < Ce^ll/L; JI<?IIl^, 



in terms of the unsealed two-point correlation matrix (2.13), for all f,g G L\ defined in (1.22). 
Proof: Windowing / and g as in ( |1.24 ), we find 



N 



N 



(3.24) {{Ll, + X)-'f,g)^, = ^ ((L^^ + A)- V., 5,)^. = E 



G\iy - x)f^{y)gj{x)dydx. 



For any p€lC, the Greens function G\ admits the Taylor expansion about Pi — Pj, 
(3.25) Ga(j/ - x) ^Gx{p,~pj) + G\{s){{y~p,) ~{x~p,)), 



for some s — s{x,y) G M. Substituting the Taylor expansion for G\ into (3.24) we observe that 

Gxip,-p,)fdy)gj{x)dydx = (®x • /)^G^ 



(3.26) 



'X- 9: 



N 

E 

« j=i ■ 

while ||G'^||l°= < Ge^. Finally, recalling the definition (1.22) of the windowed norm, we obtain 
(|3^. □ 



3.4. Residual estimates. Recalling that denotes the right-hand side of (1.2), then the residual 
J-"($) takes the form 



(3.27) 



and enjoys the following properties. 

Proposition 3.1. Let K, he a family of N -pulse configurations, as constructed in Section 2, with 
minimal pulse separation £. Then there exists r = r{£), which grows at an 0{l) exponential rate in 
£, such that for all p ^ K,, the residual satisfies the following asymptotic formula 



(3.28) 



J-i($) 

:F2m 



Ef=i (<i>r-<zr; 



hf^ + 0(e'') 
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for r — r{£) > large. Moreover, there exists C > 0, independent of e and p (z ICi such that the 
following estimate holds, 

(3.29) e (II J-i(<f>)||Li + ||Vp-^i($)|Ui) + II J-2($)||l2 + \\Vj;T2m\L- < Ce. 

Proof. We first examine J-2{'^) in the norm. Using (2.2) to ehminate the dominant terms we 
find 

(3.30) 



l-^2(<i>)|| 



L2 



< 



' I 



L2 ■ 

The first term is dominated by tail-tail interactions between the localized pulses Since these 
pulses decay at an 0(1) rate and since a22 > 1, 4'j4'k — 0(e'') for j ^ A: in any polynomially 
weighted norm where r ~ r(£) grows exponentially in pulse separation i. It follows from the L°° 
bound, (3.19), on $i that this term is 0{e^) where r can be made as large as desired. It is the 
second term which is dominant, which from the last estimate of (3.20) is 0(e). 

To examine the norm of J^i($) we use (2.8) to rewrite the first of component of J-{^) and 
apply the triangle inequality. 



(3.31) ||.Fi(<I>)||^, <£- 



|sf.i($r 



9, )0 



ai2 
3 I 



-6-1 11$^- ((Sf^l0,)"^ 



The second term is dominated by the tail-tail interaction and is 0{e^), while the first term is 0(1) 
by (3.20). The asymptotic formula (3.28) follows by identifying the leading order terms. 

The estimates on the Vp- terms follow from the bounds of ( |3.20 ) which show that ||Vp $i||l°° is 
of the same order as — qj)(f>j\\L°°. □ 



4. The Linearization and its Spectrum 
The linearization of J- about <l>(-,p) is the linear operator Lp defined by 



ail — l(j)Q!l2 



(4.1) 



a2i$^2l-i$a22 - 1 +Q22$i"$2''~ 

A direct study of these operators is complicated by the spatially varying potentials. In what follows 
we show that the singular nature of the term Lli allows us to approximate the spatially varying 
terms with finite-rank operators. In this vein we introduce the reduced linearization 



(4.2) 





° ) 




' Jll 


J12 


J21 


-^22 / 




. 






where Jn and J12 are finite rank operators 

(4.3) Jll =aiif^ . 

(4.4) J12 =012^"^ ' 
and the operators on the second row are 



(*2 -1 



ai2<j,aii-l^)^ 



(^<J,Q!l2-lfl>a 



(4.5) 
(4.6) 



N 



9^-l + «22EC'(^-^':')' 



N 



J21 =o^2i^qf^^ N 
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The differences between Lp- and are large, but the difference has zero mass in each p-window 
and hence is strongly controlled by the singular structure of L^i ■ The essential spectra of and 
Lp coincide, 

(4.7) aess{Lp) = a,,s{Lp) = B := {-e-^k^ - e"fi\k e M} . 

Although we do not pursue this issue, the point spectra of the two operators are also asymptotically 
close. 

4.1. The point spectrum. Modulo the finite- rank perturbations, the operator Ljj is lower trian- 
gular. This reduced structure affords a precise characterization of its point spectrum, up to the 
analysis of an explicit N x N matrix. So long as the pulse-pulse separation distance € > in 
JCi renders the local tail-tail interactions higher order, then the point spectrum is contr olled by 



the slowly varying component U. Estimates on the resolvent of Ln are given in Lemma 3.1 the 
following Lemma gives bounds on the resolvent of L22. 



Lemma 4.1. Fix a contour C of type (1.6) and £ > 0. Then for all (3 > Q there exists C > such 
that 

(4.8) 11(^22 -A)-Vllm^<qi/IL;,-, 

for all A on and to the right of C and for all p G /C. 

Proof. The Green's function G22 for £22 ~ A decays at an 0(1) exponential rate uniformly for A 
to the right of C since A is an 0(1) distance from the essential spectrum of Z22. The rate is also 
uniform in p since the A^-pulses decay exponentially outside 0(1) intervals. Introducing 

9 = G22 * /j we may decompose it as g = ^ gj, where gj — G22 * Xjf- From classic convolution 
estimates we have 

(4.9) 111(1 + |x -p,f).g,|Ur <C||G22|Ui|l(l + |x-p,|^)x,7IUi<qi/Li ^ 

The derivation of the bound on g' is similar. □ 

Using the bounds on the resolvents of the diagonal elements of Lp-, we establish the following 
result. 

Proposition 4.1. Fix the pulse-pulse separation distance £ > so large that the localized tail-tail 
interaction is 0{e^) for r = r{i) > 2. Then there exists > such that for all p G /C, 

(4.10) apiLp) n {5R(A) > -j/q} - (To(p) U (7s.(p), 

where (To(p) consists of N semi-simple 0{e'^) eigenvalues which are in cJp{L22), and the semi-strong 
spectrum is comprised, up to multiplicity, of the solutions to 

(4.11) ass{p) = {A|dct(/+A6.(p)) -0}, 



where the N x N matrix M\{p) is given in ( 4.20[ ). The eigenspace associated to (Tq is contained, 
up to 0(e''), within the space 

(4.12) V = span {(0, </)'i)^, • • • , (0, cj>'^)'^] . 
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Proof. The eigenvalue problem for Lp takes the form 

(4.13) -{LI, + A)*i =e-i(Jn*i + Ji2«'2), 

(4.14) (L22 - A)1'2 = - ^21*1- 

Assume that A is outside of the essential spectrum, B, and the point spectrum of ^22- Then we may 
invert L22 — A in the second equation and eliminate ^'2 from the first equation. Inverting Lf, + A, 
we arrive at the scalar equation, 

(4.15) *i - -e-i(L^i + A)-i (jii - Ji2(L22 - A)-i J21) 

Recalling (4.3)-(4.4| we rewrite the right hand side in terms of a single, finite-rank operator, 

(4.16) *i = -e~\Ll, + XyKjf ® Jr ■ *i, 
where the left and right components of the tensor product are 

(4.17) Ji 

(4.18) Jr ={an^^'''fV' - ai2$r"'(^22 ~ X)-^J2i'^>r')x- 



The eigenfunctions outside of 0-^(^22) reside inside an iV-dimensional space. To resolve (4.161 we 
act on it with (g)Jr 

(4.19) (g) • = -e-^ (g) Jr ■ {LI, + Xy\jf (g) Jr ■ *i. 
Introducing the matrix 

(4.20) Mx{p) ^ e-^ Jr ■ {LI, + , 

we see the search for semi-strong point spectrum reduces to solving the matrix equation 

(4.21) (/+AAa)(8) J^-*i =0. 



If A ^ CTp (^22) and / + J\f\ is invertible, then (gj^ -^1=0. This combined with (4.161 implies 
^'i — 0, and from ( |4.14 ), we see that \E'2 = 0. Conversely if for some v € we have {I+JV\)v — 0, 
then 

-{Li,+xr^jrv 



(L22-A)-V2i(i^i+A)-V/t/ 



(4.22) * := e^^ 

is an eigenvector of Lp. That is, X e C\i^B U (Tp{L22)j is an eigenvalue of Lp if and only if I + J\f\ 

is singular. The statement on multiplicity follows by considering perturbations of Lp which break 
any non-simple eigenvalues into collections of simple spectra and standard results on continuity of 
eigenvalues. 

To address the point spectrum of Lp arising from L22 we first recall the defining equation, (2.4|, 
for the localized pulse 0o- It is natural to introduce the linearization of (2.4 1, 

(4.23) Lo ^ dl - I + a22^o""'\ 

about (/)o. This Sturm-Liouville operator has a simple kernel, spanned by (/)q, and a single positive 
eigenvalue, Aq, corresponding to a non-zero ground state eigenfunction ^0 > 0- The remainder of 
the spectrum of Lq is contained within (—00, —21^0] for some > 0. The classical result, [l], states 
that, up to multiplicity, (Jp{L22) consists of N copies of crp(_Lo) shifted by at most O(e^), where recall 
r > 2, with eigenfunctions consisting of linear combinations of p translates of the corresponding 
eigenfunction of Lq. 
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Assume that A e (Jp{L22)\B; if the equation (4.14) is solvable, then the process follows the 
steps outlined above, and returns us to the semi-strong eigenvalue condition. The novelty lies 



in the possibility that (4.14) has a nontrivial solution with = 0, which from (4.131 requires 
that Ji2^'2 = 0. A standard implicit function argument shows that this is possible for the set of 
eigenvalues ctq := {Mk}k=i of ^22 clustered near zero (the translational eigenvalue of Lq) since the 
corresponding eigenfunctions are locally odd about each p. The corresponding eigenspace is denoted 
V. For the eigenvalues {Xo,k}k=i of -^22 clustered near the ground-state, Aq > 0, of Lq, then 



(4.24) 



N 

2k =^bkjipo{x -Pj) + 0{e''), 



and using the positivity of ipo and Xj the N x N linear system Ji2'^2k = yields only the trivial 
solution for the coefficients bkj ■ The rest of the spectrum of Lq is to the left of —2^0, and so standard 
results on the spectrum of well-separated pulses, [l], imply that the remainder of the spectrum of 
L22 lies on the real axis, to the left of —vq. □ 

We now find an explicit representation for the matrix A/a, and we show the semi-strong spectrum 
is partially characterized through the meromorphic function 

(4.25) 7^(A) ((Lq - A)- <^o""')i. , 

which is analytic except for poles at some of the eigenvalues of Lq. A similar function was introduced 
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Proposition 4.2. The matrix M\ takes the form 

E{\) 



(4.26) 

where 
(4.27) 



V A + e"/i 



0{e) 



S(A)-an'/'^-«l2«2l7^(A), 



the scaled two-point correlation function Qjq is defined in (2.12) and Q = diag((7) is the diagonal 
matrix of amplitudes. 



Proof. The reduction of the formula (4.20) for the semi-strong matrix Af\ requires the inversion of 
two second order operators. The first inversion, for L^i+X, is accomplished by ( |3.23 ). The second is 
the inversion of (^22 — A). This we achieve via the non-local eigenvalue (NLEP) machinery developed 
in 12 . For functions that are exponentially localized around the pulse positions p G JC, the NLEP 



analysis inverts ^22 — A in a windowed manner using translates of the operator Lq, 



(4.28) 



Lo,k :=a^-l + a22</'r"'(^-Pfe) 



In particular, since the potential J21 is comprised of A^-pulses localized about the positions p, the 



last estimate of ( 3.20[ ) yields 
(4.29) {L2 



N 



A)-V2i$r =«2i^gfc 
fc=i 



aii-l- 



-^(Lo,,-A)-Vr(^-m-) + (f'^). 
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where the error is in Lh^ for r > 2. Using the same estimate we may rewrite the k'th component 
of Jr, defined in ( |4.f8[ ), as 

N 

(4.30) J.fe = 5I'?r'^^("ii'^o"(^-Pfc 

fc=i 

where the function S(A) := {Lq — A) 



decays exponentially at a rate proportional to the 

A, 



distance of A to (TessiL22)- Turning to ( 4.20^ , we address the ij entry of Af) 

W'x]ij = e ^ ((^11 + ^) ^Jlii Jrj) ^2 J 

-1 {{LI, + X)-'C^,a,2a2l(^>^''-\x-p,)E{x-pJ))^, + 0{e) 



To invert L\, + A we apply (3.23) to find 



[A6 



Mi] 



X-{aii<Po'Hx-p,))Gf^®X-i. 



X- {a^2a2i<l>T-\x - p,)S{x - p,)) G^, ® x ■ ^^ + 0{e). 



Recalling the scaled, two-point correlation matrix CJat, from (2.12), and the mcromorphic function 



TZ introduced in (4.25), we may represent M\ as in (4.26). 

We are now in a position to prove the Admissibility Proposition stated in the Introduction. 



□ 



Proof of Admissibility: We denote the eigenvalues of A/a by {^J-\^\\,py\¥^i- The semi-strong 



spectrum is precisely the set of A for which ^^^\x,p] 
a single pulse, = 1, the matrix J\f\ is a scalar and 



(4.31) ^,^^\X) 
where we have introduced the constant 



-E{X)ql 



- 1 for some j — 1, 
-0(e), 



, A^. For the case of 



1 



>o Va^. 



corresponding to the amplitude of a single pulse. Recalling the operator (4.231 and the pulse 
equation (2.4 1, we find the identity 

(4.32) Lo0O = ^0 -<f>0+ + ("22 - l)0o" = («22 

which permits us to evaluate the mcromorphic function 7?. at A = 0, 



i)0r: 



(4.33) 

In particular, E{0) = 



6o" 



'L2 



a22 — 1 



^Q^-^ 7^ if 7^ 0. Moreover for |A| sufficiently far from the essential 



7^(o) = (Lo-Vr^</'o""') 

'6] 

spectrum, Tl{X) tends to zero so that E{X) — > aii^)^'-^ , while \k\\ grows, so that fi\ '{X), given m 
(4.31 ), tends to zero in this limit. As a consequence, the assumption that the one-pulse is spectrally 
stable implies there exists v e (0, ^) and & > sufficiently large such that /x^^^ is uniformly bounded 
away from —1 on the contour C^^f Moreover, since £^(0) ^ and E varies at an 0(1) rate in A 
away from its poles, there exists C > such that 

-1 



(4.34) 



1 



M^^(A) 



< C 1 
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for all A e C. 

Consider an A^-pulse and a set K. of pulse positions satisfying (1.10) for some 5 > 0. Fix a 
contour C for which (4.34) holds and consider A on and to the right of the contour. If in addition 
A is sufficiently far from the branch point — e"/i so that |A + e"/u| > ^ then 3f?(fc>,)Apj > Ce'^°' . 
Indeed, for these A the definition (2.12) and the formula (2.15) yield the estimate 



(4.35) 



Gn{X,p)Q' 



9oo '-Ny.N 



0{5 



,exp 



-6t 



and from (4.26) we see that the matrix J\[\ approximately diagonal with eigenvalues iJj^^{p,X) 



+ 0{S ^) for j = 1, • • • , N. It is here that we require S sufRciently large, independent of e, 
so that none of the eigenvalues /^j^'' attains the value —1 for A on this set. This restriction on the 
strength of the semi-strong interaction precludes the point spectrum crossing the imaginary axis 



away from the origin, thereby inducing oscillatory instabilities, as is known to happen 13 15 for 
N pulse configurations with pulses that are too close. 

On the other hand, if |A + e"/i| < e^, then E{X) is close to E{0), and is uniformly bounded away 
from zero. Moreover, one can show from an inductive proof that 



> 



N-1 



\kx\Ap, 
+ \kx\Ap, 



> CS 



N 



since A is close to the branch point, and Api > (Se~^^+"/^'. Similarly, Q is diagonal, with positive 
entries that are uniformly bounded away from zero and infinity, thus the determinant of Q^~^ is 
bounded away from zero. Moreover, the eigenvalues of Qj^Q^^^ are uniformly bounded from above, 
while we have shown that their product is uniformly bounded from below in modulus; it follows 
that each eigenvalue is uniformly bounded from below in modulus. We combine these facts into the 
estimate on the eigenvalues of N\ , 



(4.36) 



»1, 



for j = 1, . . . , iV for all A satisfying |A + e"/i| < which are on and to the right of C. Combining 
the two types of estimates we see that each A'j^'' is uniformly bounded away from —1 for all A on 
and to the right of C and for all p € JC. Moreover, the matrix / + M\ is uniformly invertible on 



these sets, in particular it satisfies (1.7), which establishes the admissibility of /C. 



We remark that since 0o and ^fp, the ground state eigenfunction of Lq, are both positive functions, 
it follows that E{X) — —oo as A — Aq from the right along the real axis, where Ao is the ground- 
state eigenvalue. In particular, if E{0) > then it must be that fj,'^^\x) attains the value —1 on 
(0, Ao) C M and hence there is a semi-strong eigenvalue on that segment. Thus the assumption that 
the one-pulse is linearly stable requires that 9 < 0. □ 



4.2. The Spectral Projection and Adjoint Eigenfunction Asymptotics. It is essential to 
our analysis that we control the projection onto the A^-dimensional eigenspacc of L associated to 
the eigenvalues ao that are algebraically close to the origin. The spectral projection takes the form 



^ (if ^^■] 

(4.37) TTffU = y ^ ' V 
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where V := {"^jjf^i and := {^-U^i are bi-orthogonal bases for the eigenspaces of L and its 



adjoint respectively. The complementary projection is given by 



TTfjU --- 



U - TTpU. 



(4.38) 

The basis elements of V satisfy the expansion (4.12); in the lemma below we develop asymptotic 
expansions for the elements of . 

Lemma 4.2. Let K, he a collection of admissible pulse configurations. Then there exists C > 0, 
independent of e > and p E JC, such that the basis elements of the adjoint eigenspace satisfy 



I* 



(4.39) 
for k = !,■■■ ,N. 

Proof: The adjoint operator is given by: 

-Lli J21 



(4.40) 

where 
(4.41) 
(4.42) 



Lt = 







L 



22 



TfhQl2-l(j,a 



4i 
jU 



For the algebraically small eigenvalues A G cto: the operator L22 is almost singular with kernel 
approximately spanned by {(f>'i^}^^i- This permits us to normalize the second component of the 
k'th basis elements as 



(4.43) 



t 

2.k 



^ -^22 '^12^1,fc 



e 142^^2242*1 



fc' 



where 7i-22 is the projection off of the small eig ensp ace of the self-adjoint operator i22- Using this 
form for j, and proceeding as in Proposition 4.2 we solve for the first component of the adjoint 
eigenvector 



(4.44) 



= [! + e-^irr Jr(/ + <)-! ® jl-]L^{J,,cj,'„ 



where J I = J/ = ^, 

(4.45) Jt := (an1>r'*i""' " ai2J2ii2"2'^22<l'2 X, 
and the N x N matrix A/jj is given by 

(4.46) =e-^®j}- L^^Jr. 



It remains to bound the second term in (4.43 
address the latter, where employing the bound ( 



\LiiJ2l'P'k 



kWw' 



< C 



(.1 



and the whole right-hand side of (4.44 1. We first 
3.4 1, we find 

X- J2l(l>'k\ - 



\J2l(l>'k 



(4.47) 
with 

(4.48) J2i0'fc = a2igr^"Vr0fe + 0{e^). 

Due to even-odd parity, this term has algebraically small mass and we deduce that 

(4.49) J2i</'fc| =0(e'-). 
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On the other hand, J2i4>'k decays exponentially away from x — pk and hence ||^2i</'fc||Lj ^ — 0(1). 
As a consequence of these two facts we find that 

(4.50) \\L^,^J2icl>',\\^i.i < Ce\ 

However, is uniformly bounded in since the small eigenspace of L22 is projected away by 7:22- 



Using (3.11 we deduce that 
(4.51) 



e-'L^^Jf\\^^.^ <C\\jr\\L^ <C. 



Finally, the matrix / + = [I + A/q)^ is boundedly invertible since /C is admissible. Taking the 



' norm of (4.44), we conclude that 



(4.52) 



The remainder of (4.39) follows by applying the estimates above to (4.43) 



□ 



5. Resolvent and Semi-Group Estimates 

In this section we generate resolvent and semi-group estimates for the reduced operator Lp>, where 
we have chosen an admissible class /C of semi-strong A^-pulses with associated contour C C C, as 



defined in ( 1.6 ) 



Lemma 5.1. For A G C, any fixed F — {fi,f2)'^ G ^^(R-); o,nd \ E C, denoting the action of the 
resolvent of L on F by {L ~ X)^^F = (91,92)'^, then 



(5.1) 

and 

(5.2) 



91 



= (L^i + {e-\jf {I + NxT^ (» Jr ■ {Ll^ + - /) KF. 



52 = ^22 - A (/2-^2l5'l), 



where Ji, Jr, and J\f\ are defined in (4.171, (4.18), and (4.20) respectively, and 
(5.3) KF = fi-e-\h2{L2 



A)-V2. 



Proof. We recall the form (4.2 1 of L and write the resolvent problem as (L— A)(5i, 52) = (/i, /2) 



For A S C the operator L22 — A is invertible. Solving for 52 and following the derivation of (4.161, 
we rewrite the first equation as, 

(5.4) (L^i + A).gi + e-\jf ® Jr ■ gi = -KF. 

We invert + A and then project with the finite-rank operator ®Jr, to obtain the matrix system 

(5.5) ®Jr-gi+ ® Jr ■ (ill + ^r^Jf ® Jr ■ gi ^ ~ ® Jr ■ (^n + X^^KF. 



Recalling the matrix A/a from (4.20), we may re-write this expression as 
(5.6) {I + Mx)®Jr-gi ^-(E)Jr-{Lli+Xr^KF. 



Since A ^ ass we may invert / + A/a to solve for (gjj^ • gi. Substituting this expression into (5.5) 
and isolating gi, we estabhsh the closed form expression (5.1 1. □ 
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5.1. Resolvent Estimates. Let K be an admissible set of N pulse configurations and let TTp- be 
the projection onto the small-eigenvalue eigenspace, (To, of Ljj, see (4.37), and let tt^ denote the 
complementary projection. We recall the norm || • \\x introduced in (1.251, and the subspace X^; 
corresponding to the range of t^jj. 

Proposition 5.1. There exists C > such that for all X in C, F e X^;, and p & JC, we have the 
following resolvent estimates for Lp, 

e 
e 



(5.7) 



(5.8) 



\\iL-X)-'F\\x < C, 
\\{L-X)-^F\\x < C, 



ell/ilU, + II/2IIL1 . 



A\h\\wy + \\f2\\L^ . 



mx) \ \kxv 

If in addition the coarse-grained projection of fi is small, then we have the enhanced residual 
estimate 

(5.9) \\{L-X)-^F\\x <C 

Proof: Using the notation 
(5.10) 



iL-X)-'F = 



e|®X-/i|+e1l/i||L;.+ ll/2||Li . 



.91 

92 



we apply the ' norm to gi as represented in (5.1 1 and use the estimate (3.1 ), 



(5.11) \\g,\\^,., <C^^\\jf{I + Mxr'®Jr-{Ll, + X)-'KFh^+C^^^ 



-\\KF\\ 



Since IC is admissible, we have the estimate (1.71 on / + Af\ for A £ C. Together with the uniform 
bound on || J;||li = we have 

(5.12) \\JF{I + Af,)-'h^<c(^l + j^^ \ 

Contained within (E)Jr is (L22 — A)~^ which is uniformly invertiblc from to for A G C since 
F S Xjj. Together with ( |3.2[ ), these observations afford the estimates 

2 

(5.13) I ® J, • (L^i + Xr^KF\ < II J.llLjKiii + X)-^KF\\l^ < C^\\KF\\li, 
Combining ( |5.12[ ), ( |5.13[ ), and ( |5.11[ ) we find 

(5.14) \\9i\\w}-^<C^\\KF\U.. 



Estimating the right hand side, we have that 

(5.15) \\KF\\li < II/iIIli + e-^ll Ji2(i22 - A)-V2||li. 
Furthermore, 

(5.16) II Jl2(^22 - A)-V2||li < II Ji2|1l2|1(L22 " A)-V2||l^ < CUMU^, 

which leads to the bound 

(5.17) IMw^:^ < (ell MU. + 11/2IUO ■ 
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Taking the H}^ - norm of both sides of (5.2) and applying (4.8), we obtain the bound 



(5.18) 



115211^1 11(^22 - A) ^{f2 - J2igi)\\H^ ^ < (^(11/21^1 ||^2lffl||Li J, 

< ciWhh^ ^ + \\9i\\wy)- 



Combining (5.17) and (5.18) we obtain (5.7). To obtain (5.8), we take the W^' norm both sides 
of (5.1 ) and then split the right-hand side into two parts, 

(5.19) <C^^\\jfiI + M^r' ® Jr ■ (LI, + X)-'KF\\L^ + \\{Ll, + X)-'KF\\^..., 



where we applied ( |3.l[ ) to obtain the second line of the inequality. We address the first factor of 
the first term on the right-hand side, 



(5.20) 



Jr ■ (LI, + X)-'KF\ < \\Jr\\L^\\{Ll, + X)-'KF\\l^ < CWiLl, + X)-'KF\\^i,i. 



Recalling both (5.12) and (5.18), while usi ng (|3.1[ ) to bound the KF term on the right-hand side 
of (5.20), we find (5.8). For the estimate (5.9), the bound on the /2 component is as before; so 
without loss of generality we consider the case F = (/i,0)^. Taking the W^'^ norm of (|5.1|), we 



have 



(5.21) llgill^r.r < C (e-'Wm, + X)-'jf\\^u^ |® J, • (L?i + A)-Vi| + IK^n + X)-'h\\^^,^) . 
Using (3.5) and the uniform bound on Jr, we obtain 

(5.22) I ® J, • (LI, + A)-Vi| < II Jr||Li||(i?i + A)-Vi||l^ < C (-^\^x-fi\ + e'll/iL; 



From (3.1 ), we have the bound 
(5.23) 



\m,+x)-'jf\\^y<c 



\j, 



1 < c 



Finally, applying (3.4) to the remaining term, 

/ £2 



(5.24) 



\\{L\, + \)-'h\\^.y<C 



ll/l 



Combining these estimates, we have (5.9). 



□ 



5.2. Semi-group estimates. Let K be an admissible set of A^-pulse positions with respect to the 
contour C. For fixed /C, the bounds on the point and essential spectrum of Lp, together with the 
resolvent estimates of Proposition 5.1 show that the operator is sectorial, and from classical result, 
, we can generate the semi-group, Sp{t), from the Laplace transform of the resolvent of Lp. 



33 



Except for the N point spectrum eigenvalues, ctq, the spectrum of L lies entirely on the left-hand 
side of C and the semi-group associated to Lp is given by the contour integral 

1 



(5.25) 



Sp{t)F 



2ni 



e^\X- Lj;)'^FdX, 



for all F e Xff. The following estimates hold on the semi-group. 
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Proposition 5.2. For any to > 0, there exists C > such that for all p £ /C^^j^, F e Xp, and 
t>to the semi-group satisfies 



(e|l/l|lL,+|l/2|l 



T.P , 



wi-i + ./2 Li . 

S T.P 



< Ce- 



InelllFl 



(5.26) \\S{t)F\\x <Ce'' 

(5.27) \\S{t)F\\x <Ce-^"''*(|lne|||/i| 
// in addition the coarse-grained projection of fi is small, then we have the improved estimate 

(5.28) \\S{t)F\\x < Ce'^^"' (e| ® X ' M + e'll/ilL; . + IIMIl^ ) ■ 



Proof. Applying each of the resolvent estimates (5.7)-(5.9) from Proposition 5.1 to (5.25) and using 
the integral bounds (5.29) we find ( 5.26 )-( 5.28) respectively. □ 

The following lemma affords estimates of the integral of the resolvent over the contour C. 



Lemma 5.2. Fix pi,p2 > 0, to > and let the contour C be as in (1.6), then there exists C > 
such that for all t > to 

1 a pi + P2 < 2, 



(5.29) 



mkx)\'"\kx\^' 



|lne| if pi+p2 = 2, 

£-f(pi+P2-2) if p^+p2>2. 



Proof. Estimates on the angled parts C/± are straightforward because of the exponential decay in 



The concern is for the vertical part, where for A = — e"t/ + is, we have 



(5.30) 



for some constants C± > 0. Applying these estimates (5.29) we have 
(5.31) 



dX < Ce-'^^'e-fPi+P^) 



ds 



(s2 + e2°iy2) 



P1+P2 



For pi + P2 < 2 the integral is uniformly bounded, even if /i = 0. For pi + p2 > 2 we rescale s by 
e" and bound the resulting integral. □ 

Remark 5.1. The semi-group estimates are derived for t > to- Short time estimates can be 
derived which scale like t~^^'^ as t ^ 0. This integrable singularity has no impact on the analysis 
which follows, and we omit it. 

6. Nonlinear Adiabatic Stability by Renormalization Group 

Our primary result is that, in an || • ||x neighborhood of each admissible A^-pulse manifold /C, we 
may decompose solutions U = {U, V)"^ of (1.2) as 

(6.1) U ^<^p-^W*{x,t), 

where the iV-pulse configuration p = p{t) is dynamic in time and the remainder W* = (W"i , )^ 
can be uniformly controlled in the || • ||x norm. In particular, we prove the adiabatic stability result, 
(1.14), and derive the limiting pulse dynamics, (1.15). The pulse profiles, $, are only approximate 
equilibria. Inserting the decomposition (6.1) into (1.2) and Taylor expanding J- about (^p yields 

(6.2) W* + Vp-$ • J'($) + LpW* ^ X{W*), 
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where the hnear operator Lj^ was mtroduced in (4.1 ) and the nonhnearity X — {Xi, ^2)^ is 

(6.3) X{W*):^F{<^ + W*)-F{H-LpW*. 

The renormahzation group procedure starts by freezing the reference pulse configuration p = po 
in Xpg, where po is as constructed in Proposition 6.1. The semi-group estimates of Proposition 5.2 
require a hnear operator with a frozen pulse configuration, so that we replace Lp not merely with 
its reduced linearization, but with its reduced linearization at a frozen pulse configuration, Lp^. 
Moreover, we separate the remainder W* into two parts, W* = W+^c{-,P,Po), where the correction 
term, $c, which serves to adjust the shape of the pulse profile, depends upon both the fast pulse 
configuration, p — p{t), and the frozen reference po. More specifically, we have the decomposition 



(6.4) 



U = ^p+^cip,Po) + W, 



where the correction term, <i>c = {^0,1,^0,2) is chosen to cancel the parts of the residual, 
which do not contribute to the pulse dynamics, 



(6.5) 



$c(-;P,Po) -LJnp^T{<^i-,p)) 



Here Lp^ is the reduced operator defined in (4.2), and is the orthogonal spectral projection 
off of its small eigenvalue eigenspace, (To(po)i defined in (4.38). The complementary part of the 
remainder, W = {Wi, W2)'^ , incorporates errors that arise from the transient nature of the flow - 
parts which may not necessarily be slaved to the pulse configurations. 

The following proposition constructs a base point po about which the local coordinate system is 
developed. 

Proposition 6.1. Let IC be an admissible class of N -pulses. Then there exist 5 sufficiently small, 
Ml > 0, and a smooth function J : X ^ JC such that for all Uq — {Uq, Vq)'^ satisfying 

(6.6) \\Uo - {Hp.) + $c(p*,P*)) llx < <5, 

where := Uq — {^{p*) + ^c{p*,P*)), for some E JCs, then p = p^, + J'{W^,) satisfies 

(6.7) Wo Uo - {Hp) + <^c{p,p)) e Xp, 

for Xp defined in ( |l.ll[ ). Moreover, if W^ :~ Uo ~ (<E>(p*) + $c(p*,p)) G Xp for some p E ICg, then 

(6.8) \p-p*\<Mi\\W4x\p*~p\- 

Proof We may write Uq = $(p*) + ^c{p*,P*) + W^ and Uq = $(p) + ^c{p,p) + Wo, and hence 

(6.9) Wo^w, + Hp,) + ^c{p„p,) - {Hp) + '^Mp)) ■ 



Since by (6.5), 7rp^c{p,P) = 0, the equation (6.7) is equivalent to solving the nonlinear system 



(6.10) 



A,{p,W,;p,) (w,+Hp,) + H{p*,P*)-Hp),'^]{p)) 



L2 



Writing A = (Ai, • • • , An)'^ , it is clear that A(p*, 0) 0. In hght of ( |3.19[ ) and ([7^ the function 
A is smooth in p, p, and W, . We examine the following gradient 



3.11) 



We evaluate the right-hand side using (3.19|), (|3.20|), and (|4.39|), and obtain 
(6.12) ypA\p=p„w,=o = 



-2021/(022-1) _|_ 0{e), 
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where Q is the diagonal matrix of q amphtudes, whose determinate is bounded away from zero 



for p G /C. The existence of J verifying (6.7) fohows from the impHcit function theorem. If, in 



addition, we have W^, := Uq — {^{p*) + <i>c(p*,p)) G Xp, then we may re-write (6.101 as 
(6.13) 



where the last equality follows from the assumption W^. € Xp. We use (6.12) to Taylor expand the 
left-hand side, obtaining the leading order relation 

^ ^ 1 



(6.14) 



However, we readily see from Lemma 4.2 that |^| < C||M^*||x|p — p|, which establishes (6.8|. □ 



6.1. Projected equations. We emphasize that the pulse configuration p — p{t) in $ = $(-,p(i)) 
is not frozen, but evolves freely with t. The separation of the pulse configuration into fast and slow 
(frozen) variables requires the introduction of the secular operator, AL := Lp— Lp^, which induces 
a temporally growing forcing term as p — p{t) evolves away from po- Eventually the secular term. 



p{t) — Po, grows too large, and the equations are re-projected, using Proposition 6.1 to update the 



frozen base-point po to pi and Wq € Xp^ to W\ G Xp^ . In the sequel, we show that the long-time 
evolution of both pulse configuration and remainder is well described by the collection of initial 
data {p„}^^o and {W^„}Jf=o- 



As the first step in this process, we rewrite (6.2) as 



(6.15) 
(6.16) 



Wt + (Vp-$ + Vp-^c) ■P = ^p-o^($) + Lp;^W -f AL ($e + W^) + A'($e + W), 

iy(a;,0) =1^0, 



where pq and G Xp^^ are as determined in Proposition 6.1 Moreover, following classical mod- 
ulation theory, we choose the evolution oi p — p{t) to enforce the condition W G Xp^. To ex- 



tract the evolution of the pulse configuration, we project (6.151 onto the small eigenvalue adjoint 
eigenspace, {'^]{-;Po)}jLi, whose orthogonal complement coincides with the kernel of npg. Since 
0, we attain the fast pulse evolution 

Mp = g{p;po,W), 
p{to) ^Po, 
where the N x N matrix M takes the form 

dpi Upi - / ^2 

and the forcing term g G is given by 

(6.19) g, (^($(p)) + Ai ($, + M^) + A'($, + T4^), *](po))^^ = (j-($ + M^*), *](po))^^. 

The evolution of the remainder W is obtained by applying the complementary spectral projection, 
TTp^, to (6.15), which yields 



(6.17) 



(6.18) 



IS'*'"* 



(6.20) 



Wt = R + Lp^,W + np„{AL{^, 
W{x,0)^Wo, 
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where we have introduced the temporal component of the residual, 
(6.21) R{f;pa) (Vp-$ + Vp-$c) • ^■ 

The central goal of the renormalization approach is to control the growth of the remainder W 
and the secularity, \p{t) — po\, in systems (6.17) and ( 6.20[ ) in terms of the quantities 



(6.22) 
(6.23) 



Tw{t) '■— sup 



^e°i/(s-to)| 



to<s<t 



W{.s)\\ 



X, 



Tp{t) := sup \p{s) - Pol- 

to<s<t 



6.2. Bounds on Pulse Dynamics. The following lemma verifies that the matrix M is uniformly 
invertible. 



Lemma 6.1. The matrix M defined in (6.18) has the following asymptotic form 



M = - 



(6.24) 

where Q is the diagonal matrix of amplitudes q. 



Proof. From (4.39) we see that 5'' is dominated by its second component, which satisfies 4* 



0^ + 0(ei+"/2). Similarly, from ( |3.20[ ) we see that ^ = -(?!)^ + 0{e^+°'/'^). From we have 

the bound HVp-^dlx ^ Ce^^"/^, which shows the contribution from this term is higher order. 

Recalling the scaling (2.3) yields (6.24). □ 

The lemma below establishes an upper bound on the rate of pulse motion. 



Lemma 6.2. The pulse evolution, given by (6.17) satisfies the following bound 
(6.25) bi <c[e+{e + r,-(t))(|| W^(.)||x + e^""/^) + 11^^)111 + e^^) . 



Proof. From (6.24) the matrix M is boundedly invertible. We apply — L°° estimates to the first 
component of g and — estimates to the second component, 

(6.26) I < C (ll J-i($ + W*)\\l^ ll^iilU- + ll-^2($ + W^*)|U2||vI/t 



for T{<S> + W*) = (J'i($ + W^*), J"2($ + W*)y defined in (|6l9|. From (|4^ we have ||^'].ilU- < 
Ce^ while H^j^slli^ = 0(1). Combining these estimates with ( |3.29[ ), ( frsj ), ^L6^ , ([T^Sf, and the 
bound \\W*\\x < C{\\W\\x + e^-^/^) afforded by ^ yields ( [6^ . □ 

6.3. Decay of the Remainder. The following proposition establishes uniform estimates on the 
decay of over the duration, At, that the linearized operator L ~ Lp^ is fixed at pq. It is 

convenient to introduce the decay factor 

(6.27) e-'''''^\ 

which bounds the action of the semi-group Sp^ (At) on Xp^ . 

Proposition 6.2. Fix ep > sufficiently small and let the normal hyperbolicity condition, < a < 
i, hold. There exist constants Ci, C2 > such that for all < e < and all initial data Wq G Xp^ 
satisfying 



(6.28) 



In el 
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then the solution W to \6.20 1 satisfies 

(6.29) \\W{t)\\x < C2 (||W^(to)||x|lne|e-^°''(*-*°) + e^-") , 
for all t Cz [t(),to + At] for any At for which the decay factor meets 

(6.30) l3{At) > Cie^-i". 

In particular, we may choose At so that at ti := to + At we have 

(6.31) \\Wit,)\\x < C2(||W^(to)|U|lne|ei-3"/2 + ei-") , 

(6.32) \p{ti)~Po\ < C2\\ne\e^-". 

Proof. Applying the variations of constants formula to ( 6.20[ ), we have 

(6.33) Wix,t) = Sp-,{t - to)Wo + J^^ 5p-o(i - s) (i? + (AL($, + W) + <¥($, + iy)))ds, 



where Sp^ is the semi-group generated by Lpg. Applying (3.191, (3.201, and (7.11 to (6.211, we find 
(6.34) e"/2||i?i|Ui + ||i?2|Li _ <C\f}\. 



We take the X-norm of (6.33), and apply the semi-group estimate (5.261 to R and (5.27) to W[ 



while using (7.14)-(7.16) and (7.1), on the remaining terms. Powers of ||M^||x appear from several 



terms, but the dominant contribution is from the temporal residual, \p\, via (6.25), yielding 
(6.35) limOllx <c(^e-^°'^(*-*»)|lne|||iy(io)||x+^'e-^°''(*-^)|p'|rfs^ . 



We evaluate (6.35) at t = t', multiply by e^""^*^* use ( |6.25 ) to controlp and bound ||W(s)||x < 
g-e°iy(s-io)y^(^^^ vahd for t > s. Taking the sup over t' e (^o, ^), we obtain 



Tw{t) < C\\ne\Twito)+C e^^^f^"*") (e + e^-") + 

Jto L 

(6.36) (e + Tp-(t)) (Twit) + ei-"/2e^°K^-to)^ + T^{t)e-'''''('-*°^' 



ds. 



Evaluating the integrals on the right-hand side, recalling the decay factor /3 from (6.27), and keeping 
dominant terms, yields 

(6.37) Twit) <c(^\ \ne\Twito) + ^ + + (Twit)At + "^^^ + ^""^^.(i) 

where At :~ t — tg is the length of the renormalization interval. The following lemma bounds Tp. 

Lemma 6.3. Fix /C, then for a < 2 there exists a constant Ci > sufficiently small, but indepen- 
dent of e, and C2 > such that for all t > tQ for which 



(6.38) 

then 

(6.39) 



Twit) < Cie" and At < Cie"/^-!^ 

T>(t) < C2 (Ai(e + e^-") + e^-^Twit) + £""^^(0) • 
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Proof. From the definition (6.23) of Tp- we have the bound 

i-to+At 

(6.40) Tp{t)< / Wis)\ds. 



Using the estimate ( |6.25[ ), the bound |lM^(s)||x < CTwit)e-^°'''^''~*°'> vahd for to < s < t, and 
integrating in s yields 

(6.41) Tp{t) < C (Ai(e + e^"") + e'-"Tw{t) + T,j{t){e-"Tw{t) + Ate^-^/^) + e-"T^{t)') . 



Subject to the constraints (6.38), for Ci sufficiently small, the term on the right-hand side of 
(6.41) may be absorbed into the left-hand term. Adjusting the constant we obtain (6.39). □ 



We continue the estimation of T^r required to establish Proposition |6.2| Inserting the bound 



(6.39) of Lemma 6.3 into (6.371 and collecting powers of Tw we find 

/ / fl-" + (At) (e2-3a/2 ^ g3-5a/2^ 



Tw{t)<C\ \\ne\Twit, 



0) 



/3 



+ 



Twit) ( ^^4^ + eAt + (Atfie + e^-") ) + 



(3 



(6.42) 



1-5q/2 

"2 '^-^ I + ei-"At ) + Ate-"T^{t) 



/3 



For this estimate to be meaningful we require that the coefficient of Tw on the right-hand side be 
strictly less than one. This generates two conditions, the first requires 

2-5a/2 1/5 

(6.43) 13 > — - — , equivalently At < e""- (2 - -a)| lne| - InC 

and since the renormalization interval At must be at least e^" in order to insure that /? <C 1, we 
impose the bounds 

(6.44) e"" <At< e-"^ (^(2 - ^a)| In e| - In C 

In light of these bounds, the second condition for the absorption of the T^, term becomes 

(6.45) {Atf (e + e^-") < i equivalently 



1 

a < - 
2 



which sets the ultimate limit on a. Subject to these conditions we absorb the linear term and merge 
the cubic term into the larger quadratic term. Enforcing the bounds on At we obtain the key 
estimate 

(6.46) Twit) <c(^\ \ne\Twito) + ^ + ^lit) (^e"" + '-^^ 

This inequality has the form of a quadratic equation, h(r) = where r Tw{t) and 

an ^2 



(6.47) 



hir) ■.^C[\\ne\Twito) + 



-r + C{e 



l-5a/2 



13 
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This equation has two positive roots < ri = O{ao) ^ r2 if 0002 ^ 1. Indeed, taking Ci in (6.28) 
sufficiently small, independent of e, these are precisely the conditions imposed by (6.28) and (6.301, 
under which we have the bound 



(6. 



Tw{t)<ri<C[ \ \ne\Twito) 



Recalling the definition (6.22) of Tw we obtain (6.29). Choosing the minimal value of /3, equiv- 
alently the maximal value of At, permitted in (6.301 we find the first bound of (6.311, while the 
second follows from (6.28). The estimate (6.32) is then a consequence of (6.39) and the bound 
dels! on Tw □ 

6.4. The Renormalization Group iteration. At the conclusion of the first renormalization 
interval (to, ^i), where ti :— to + At, we have a pulse configuration p{ti) and a remainder W{-, ti) G 
Xp-Q whose X norm is smaller than that of W{-,tQ), so long as a < i. We are in a position to 
iterate the renormalization process until p hits the boundary of IC. That is given Pn-i € ICq, 
tn :— tn-1 + At„_i, W{-,tn) G -'^p„_i Satisfying (6.28), and p'(t„_i) satisfying (6.32) then defining 



from Proposition 



6.1 



?7(-,t„) = $(-,p(t„)) + $c(-;p(i«),Pn-i) + W^(-,i„), 
we may construct p„ p{tn) + J{W{tn)) and VF„(-) G Xp^ such that 

[/(•,t„) = <i>{-,p„)+<Pc{-,Pn,Pn)+Wn{-), 



where <i>c is defined according to (6.5). Moreover, since W{tn) € Xp^ -^ we may apply (6.8) to 
estimate the jump function JJ, 



(6.49) 



\Pn-pitn)\ < Mi\\Wit„)\\x\Pn-Pn~l\, 



which, from (6.31 ), is negligible compared to evolution of p{t) over t e (t„_i, i„). The jump in the 
remainder upon the renormalization satisfies 



(6.50) 



\\Wn - W{t„)\\x = \\{^{Pn) - $(p(tn))) + {<i'c{Pn,Pn) ~ <Pc{Pn~l , p{tn)))\ 



X 



The dominant contribution in the first term comes from the $1 component, which from (3.19) has 
a sensitivity to variation in p on the order of e^". Using (7.1 ) on the second term yields the bound 

(6.51) ||T^„-l^(i„)||x <c(e-"|p„-p(i„)H-ei-"/2|p„-p„_i|) . 



The estimates ( |6.49[ ) and ( |6.32| yield 
(6.52) \\Wn-Wit„)\\x <C 



lne|ei-2o||^^||^ + |ln,|e2-3a/2 



In particular, applying the triangle inequality to bound ||IVn||j(: we see for < a < | that the 
dominant contribution is from ||VF(tn)||j>s:, which from (6.31) takes the form 

(6.53) \\Wn\\x <c(\\ne\e'-'"'^^\\Wn-i\\x+e'- 
We introduce the Renormalization Group map 

(6.54) 



Pn-l 



Pn 
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which produces the initial data at time t„ from the initial data for the pulse and remainder equations, 
(6.17) and (6.201, on the interval (<n-i,^n)- From a simple, linear iteration argument applied to 
(6.53), we have the estimate 

llW^nlU < C ((l lne|ei-3"/2)" \\W4x + e^' 



(6.55) 



The estimate (1.14) follows. 



6.5. Long-time asymptotics. From (6.32) we see that Tp < C\ lne|e^~" at each point of the RG 
iteration process, while < TV- Using these bounds in (6.25), we see the dominant long-time 

contribution to the pulse evolution arises from the residual, with the next largest contribution from 
the square of H^dlx- From (6.24| we may also invert the matrix M, yielding the reduced evolution 



(6.56) 



Pj = 



q 



2q2i/(Q22 — 1) 



(^2(1'),. 



3> L'^ 



o( 



\4w\\xA\w\ 



foT j — 1, - ■ ■ ,N. Substituting for J^2i^) from ( |3.28[ ) and integrating by parts, we find 



(6.57) 



(•^2($),0;)^2 



"21 



/a22 + l^ 



"22 



Oie 



2+a\ 



Using the last estimate of ( 3.20[ ) to replace <i>i with qj, yields 
(6.58) (J-2(<i>),0;-)^, : 



CK21 ^1 

Q!22 + 1 



) +0(e2+"). 



Turning to (2.6), we see that the odd derivatives of $i at a; = pj of order 3 or higher are at most 
0(e'^+"), while the even derivatives contribute an odd component to ^'i- That is, in a neighborhood 
of a; = Pj we have 

(6.59) ^[{x) = $'i(p,) + Eom{x) + 0(e3+") , 

where E!odd(2;) is odd about x — P j an d does not contribute to the inner product. Combining 
these observations with the scaling (2.3) and inserting into (6.56) shows that the pulse velocity is 



proportional to the derivative of $i at the pulse position, 

"21 



(6.60) 



: ||a22 + l 1 
»0|li = 22+l 1 



<i>[{p,) + 0{e'-",e\\W\\x,\\Wfx) 



"22 + 1 ||(pblli2 qj 

To evaluate ^[{pj) we invert Lf^ in (2.8) in terms of its Green's function Gq, ( |2.17 ), and take the 
x-derivative of the result, and rescale according to (2.3), yielding 

N 



(6.61) 



<S>[ip,)^-e-'J2lk I 4>r{x-pk)G'^{p,-x))dx, 

1 1 



k=l 



k=l 



(/)o'"(a; -pfc)sign(x -pj)e 



The Green's function is slowly varying and (j)Q is well localized, moreover when j = k the integral 
is zero due to even-odd parity. We evaluate the integrals asymptotically, obtaining the expression 



(6.62) p = eQ-'A{p)q" + 0(e^-",e||PF|U, ||W^|| 

where Q is the diagonal matrix of the amplitudes q and the antisymmetric matrix A{p) is defined 
component-wise as (1.16). It is clear that the pulse evolution is at most 0(e), hence given initial 
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data of the form (1.12), then the projected initial pulse condition po is within 0(e"|lne|) of p 
and hence its distance to dJC is of the same order as do := d{p,dlC). Consequently the time 
to arrive at dIC is O(e~^(io). Moreover evolution of the form (1.15) generically causes the pulse 
spacing to increase in time. Indeed within the weak regime the interaction reduces to repulsive 
near-neighbor tail interaction, and the semi-strong spectra is fixed, so that the time to exit the 
admissible configuration domain, JC, is infinite. This completes the proof of Theorem |l.l[ 

7. Technical Estimates 

We prove several technical estimates used in Section 6. The first involve the correction term <I>c 
in the decomposition (6.4). 

Lemma 7.1. For each admissible family JC of N -pulse configurations there exists C > such that 
(7.1) ||$c||x + ||Vp^$c||x<Cei^"/2^ 

for all p G /C, 



Proof. From the estimates (4.39), the definition (4.37) of np, and the bounds (3.29) we calculate 



so that to the orders we are concerned with, 7rp-J^(<I>) = J^($). We take the || • ||x norm of (6.5 ) and 

gi+Q/2 obtain 



use the resolvent estimate (5.7) at A = 0, for which 3fi(fco) 
(7.2) 



l^ol 



The bounds (3.29) on the residual yield the $c estimate of \7.1\ . Taking Vp' of (6.5), and observing 
that Lp^ and TTp^^ are independent of p, yields 



(7.3) 



^^'^PoVp-F($). 



The estimate (5.7), applied at A = 0, yields the bound 

(7.4) W^p^cWx < C (ei-"/^||Vp-J-i($)|U. + e-"/2||Vp-J-2($)Li ^ 



which, when combined with (3.29) yields the Vp estimate of (7.1). 



□ 



We break AL into secular, A^L :— Lp~ Lp^, and reductive, ArL := Lp^ — Lp^^ which satisfy the 
following bounds. 

Lemma 7.2. For each admissible family, K,, of N -pulse configurations there exists C > such that 
for all V = {Vi,V2Y e Xp 

(7-5) e\\[A,LV]i\\L^ + \\[A,LV]2\\lI^^ <CTp[t)\\V\\x, 

(7.6) 

(7.7) [A.Ly]i| =0, 

(7.8) £-111^-1(^^)11^. + 11^-2(^)11^; <C||F||5 

Proof. We first examine the secular operator AgL, 

(^ifr'^p-f 



(7.9) AL, 



(j)a2l-l<j)a22 



"22 



<J)a21 (1)022-1 

^po,l^po,2 ^ 
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The $2 terms are smooth and decay exponentially away from each pulse position, in particular 
there exists a C > such that 



(7.10) 



|*p.2 - $p-o,2||//i + ll*P.2 - *p-o,2||li < C\p-po\ < CTp{t), 



where Tp is defined in (6.23). From (3.20) we have the estimates < Ce while for any fixed 

/? > 0, |j$2'&i|lL°° < C. Combining these estimates yields 

(7.11) \\[A,LV]i\\li <^mp.2-'^po,2\\L4Vl\\L- 



\^P.2~'^P0,2\\l4V2\\l-) < -Tpit)\\V\\x 



The estimate on 11 [A^i V^]2|jii is similar and (7.5) follows. For the reductive operators, A^i, the 

i,p 1 1 

difference Lp^ — L^^ is large but has small mass in each x window, which permits the application 
of (5.9). The reductive operator takes the form 



(7.12) ArL 



) 



V 



-^^ Po,l Po,2 
N 



Jl2{Po) 



"22l>p-„,i<I'p-„^2 



"22 E 00 ^i^-Pj.o) 
j = l 



where Jn and J12, defined in (4.3)-(4.4), are chosen so that 

[A,LF]i =0, 



for all V e W^'^( 

Ce~^\\V\\x appearing as the first term of (7.6 



and thus (7.7 1 is satisfied. The weighted-norm bound ||[Ariy|Ui < 
I 1 i,p 

follows from typical Holder estimates. The sec- 



ond bound of (7.6) follows from the form, (4.6), of J21 and the last estimate appearing in (3.20) 



applied to ($1 — qk)<t>k P taking the values ^22 ^ 1 and a22- 

For the nonlinearity, we consider the case ai2 — a22 = 2; the other cases are similar. The leading 
order terms in the nonlinearity take the form. 



(7.13) 



X{V) 



1<J,Q11-1 



^2ViV2 + ^"''V^ + e-^'S?'^''-^^lV{ 



-1 

where the derivatives in Vi are smooth since $1 is uniformly bounded away from zero, where the 
powers C/"^i in J- are non-smooth. For terms in which the power of $2 is positive, we use the 
bound < C which may be inferred from (|3.20[). For the terms without $2 we use 



(3.21). Together this approach yields (7.8). 

Combining the estimates above with Proposition |5 . 2| yields the following corollary. 
Corollary 7.1. For each admissible class K. there exists C > such that for all t > 
(7.14) \\S{t) {np„A,LV) \\x <CTp-{t)e-'"'''\\V\\x, 



□ 



(7.15) 
(7.16) 



S{t) {np-„ArLV)\\x <Cee' 



\V\\ 



X, 



\\S{t) {np,MV))\\x<Ce-^"'''\\Vrx. 

8. Discussion: Normal hyperbolicity and adiabatic stability 

can be motivated by 



1.1 



The normal hyperbolicity condition a G [0, |) presented in Theorem 
the following simple argument. For a fixed N pulse configuration po = pli-o) with associated lin- 
earization Lp^, we obtain decay estimates for the semi-group generated by Lp^, see Proposition 5.2 
which guarantee exponential decay over time-scales t — to ^ e~". However, over this time period 
we find that the pulse positions p — pit) experience a drift on the order of |p(to + i) ~ R) I ~ £^ 
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and the linearization about the evolving pulse configuration is time-dependent, L := L^^i^^y For 
non-self adjoint operators, time-dependence in the linear operator can act as a source of forcing 
that destabilizes the underlying equilibrium, even if for each fixed t the linearization has exponen- 
tially stable semi-groups. In general, the problem of characterizing the semi-group produced by a 



time-dependent linearization is nontrivial, 37 38 . However, if the time-dependence of the linear 



operators is sufficiently slow compared to the exponential decay rate of the semi-group associated to 
each fixed-time operator, one would generically imagine that the time-dependent operator inherits 
the exponential decay of its frozen-time constituents. 

A characterization of normal hyperbolicity for a fiow in a neighborhood of a manifold should 
specify how slow is "slow-enough." We introduce a characterization via the "secularity" in the 
linearization, AL := ip-jj) — Lp-^ . Introducing the spectral projection, TTp^^, onto the tangent plane of 
the manifold of semi-strong iV-pulses at pulse configuration po , the forcing induced by the secular 
term is characterized by the quantity ||(/ — T:p„)/S.L\\ where || • jj denotes an appropriate operator 
norm. Assuming the generic scalings 

11(7 - ^pJALll ^ II ALU ^ \p{to + t)-p^\^ et, 

then over the time-scale, t — to^ e "7 required to obtain exponential decay, the contribution of the 
secularity scales like 

/ 11(7 -7rp-JAi:(s)||e-''^°(^"°-'*)ds- 6^-2" / se-'^^-'Us = 0{e^-'^°'). 
Jo Jo 

This suggests the normal- hyperbolicity constraint a < ^. In the study here-in, this constraint is 
manifest in equation (6.42), which requires (Ai)^e ^ 1, where At ~ e~" is the time required to 



obtain 0(1) decay in the semi-group. This discussion motivates the normal hyperbolicity conjecture: 
"The square of the linear decay time multiplied by the pulse velocity must be sufficiently smaller 



than one." For the system (1.2), this is the genesis of the normal hyperbolicity condition, a < ^. 

What structure remains for decay rates slower than the normal hyperbolicity threshold? Our 
intuition is that the RG iteration, which effectively produces linear semi-group estimates for the 
time-dependent linearization, Lp(t) , will fail, and the separation of the spectral spaces into dynamic 
(the translational eigenvalues) and slaved (the rest of the spectrum) modes breaks down. It may 
be possible to recover the required linear estimates if the nose of the essential spectrum, which 
extends to — e"/i, is incorporated into the dynamic elements of the decomposition. In particular, 
resonance poles, the zeros of det(/ + Af\) for A in the essential spectrum (see Proposition 4.1), 
must be accounted for. It is unclear if this can be achieved in a finite-dimensional setting, or if the 
underlying dynamics would be truly infinite dimensional. The situation is evocative of the spectral 



gaps required to produce inertial manifolds for dissipative PDEs, see 23 . It is precisely the spectral 
gap which is closing here. 
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